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ABSTRACT 

Aims. We study non-axisymmetric oscillations of a straight magnetic tube with an elliptic cross- 
section and density varying along the tube. 

Methods. The governing equations for kink and fluting modes in the thin tube approximation are 
derived. We found that there are two kink modes, polarised along the large and small axes of 
the elliptic cross-section. We have shown that the ratio of frequencies of the first overtone and 
fundamental harmonic is the same for both kink modes and independent of the ratio of the ellipse 
axes. 

Results. On the basis of this result we concluded that the estimates of the atmospheric scale 
height obtained using simultaneous observations of the fundamental harmonic and first overtone 
of the coronal loop kink oscillations are independent of the ellipticity of the loop cross-section. 
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1. Introduction 

k>" ■ The solar atmosphere is a highly dynamic and structured plasma that is able to support a wide 

■ variety of magneto-acoustic waves and oscillations. Each layer of the solar atmosphere, from the 

■ photosphere to the corona, is magnetically connected to the others via the all pervading magnetic 

_ _ i 

field. The omnipresence of the waves throughout the atmosphere is becoming well documented as 
new and exciting techniques are being de veloped to help observe and study the waves (see, e.g. 
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Baneriee et al.ll2007 



Tomczvk et al.ll2007 ) . 



After transverse coronal loop oscillations were first observed by TRACE in 1998 



(lAschwanden et alJll999t iNakariakov et alJI 19991) . the phenomenon became one of the hot top 



ics within solar physics. In the first theoretical interpretation of these oscillations, a coronal loop 
was modelled as a straight magnetic cylinder with the density constant inside and outside. Since 
then, a number of more complicated and realistic models have bee n considered. For a recent rev iew 



on the theory of transverse oscillations of a coronal loop see, e.g jRuderman & Erdelvil (120091) . 
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Fig. 1: Sketch showing the elliptical coordinate system used to describe the loop cross-section. 
The open and closed curves show the s and <p coordinate lines respectively. The thick closed curve 
shows the tube boundary. 



Although the transverse coronal loop oscillations are interesting on their own, their main impor 



Nakariakov & Of man 



tance i s related to the fact that they are a powerful tool of coronal seismology. 
(120011) demonstrated this by using the observations of transver se coronal loop oscill ations to esti 



Andries et al 



d2005l) suggested to 



mate the magnitude of the magnetic field in the corona, while 
use these observations to estimate the atmospheric scale height in the corona. 

In this paper we continue to study the transverse oscillations of coronal loops. Coronal loops 
with elliptical cross-se ctions and a constant density profile h ave been studied previously in both 
cold ( RudermanlboCB ) and finite-/? ( Erdelvi & Morton 2009 ) plasmas. Now, we consider oscilla- 
tions of loops with the density varying along the loop and a constant elliptic cross-section. The 
paper is organized as follows. In the next section we formulate the problem. In Sect. [3] we de- 
rive the governing equations for non-axisymmetric oscillations of a coronal loop with an elliptic 
cross-section in the thin tube approximation. In Sect. [4] we study the implication of our analysis on 
coronal seismology. Section[5]contains the summary of the obtained results and our conclusions. 



2. Problem formulation 

We model a coronal loop as a straight magnetic tube with an elliptical cross-section. The cold 
plasma approximation is used. The density varies along the tube, while the cross-section remains 
constant. In Cartesian coordinates x, y, z the loop axis coincides with the z-axis. The equilibrium 
magnetic field is given by B — Bz, where B is constant and z is the unit vector in the z-direction. 
The plasma motion is governed by the linearised ideal MHD equations, 

H = J_ (V x b) x B, (1) 

at 1 nop 

^Vx(fxfi). (2) 

Here £ is the plasma displacement, b the magnetic field perturbation, p(z) the equilibrium density, 
and //() the magnetic permeability of free space; p(z) = p;(z) inside the tube and p(z) = p e {z) outside 
the tube. 
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Let us introduce the elliptic coordinates s and <p in the xy-plane (see Fig. 0. The Cartesian 
coordinates are expressed in terms of elliptic coordinates as 

x - <x cosh s cos ip, y = cr sinh s sin <p, (3) 

where cr is a quantity with the dimension of length, s varies from to oo, and <p from -it to n. In the 
elliptic coordinates the equation of the tube boundary is s = so. Then the large and small half-axes 
of the tube elliptic cross-section are in the x and y-direction, and they are given by 

a = crcoshso, = crsinh so- (4) 

At the tube boundary the normal component of the displacement, g s , and the magnetic pressure 
perturbation, P = b ■ B/po, has to be continuous, 

161=0, [P]=0 at s = s , (5) 

where [/] indicates the jumps of function / across the boundary defined as 

[/] =lim[/(* + £)-/(*-£)]. (6) 

The magnetic field lines at the loop foot points are frozen in the dense photospheric plasma, so 

that 

£ = at z = +L/2, (7) 

where L is the loop length. 

It follows from Eq. (0 that the points with the elliptical coordinates s — 0, tp = (po, and 5 = 0, 
<p = -<po are the same point in the xy-plane. This implies that P and £ s have to satisfy the boundary 
conditions 

P(0, <p) = P(0, -cp), ^(O, ip) = -6(0, -ip). (8) 

Equations ((TJ and (|2]i together with the boundary conditions (0, (0 and (0 will be used in the 
next section to derive the governing equations for non-axisymmetric oscillations in the thin tube 
approximation. 



3. Derivation of governing equations 



Dvmova & Rudermanl (120051) to derive the 



The analysis in this section is similar to one used by 
governing equation for a thin tube with a circular tube cross-section. We begin by noting that, in 
accordance with Eq. 0, 6 = 0. The system of Eqs. (0 and (0 can then be transformed to 



dt 2 p /j.qp dz 

oz 

P = -pv 2 A V-£, 



(9) 

(10) 
(11) 



where v& is the Alfven speed defined by v\ = B 2 /pop, and the operator Vj_ and component of the 
magnetic field perturbation perpendicular to the z-axis are given by 

. d 



'dz 



b x =b -b -2. 



(12) 
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Eliminating b± from Eqs. (|9]l yields 

Taking the divergence of this equation and using Eq. ( TTTb we arrive at the equation for P, 



< 14) 



Now we use the thin tube approximation. To do this we note that the characteristic spatial scale in 
the z-direction is L, and the characteristic time of the problem is L/Va, where v& is a typical value 
of Alfven speed. In what follows we only consider the perturbations that decay at the distance of 
a few a from the tube. Then the characteristic spatial scale in the x and y-direction is a. It follows 
from this analysis that the ratio of the left-hand side of Eq. (TT~4T > to its right-hand side is of the 
order of {a/L) 2 <k 1, so that we can neglect the left -hand side. Then, using the expression for 
in the elliptical coordinates (e.g. Korn & Korrll961 ). we obtain the equation for P in the thin tube 
approximation, 

— + —=0 (15) 

ds 2 d<p 2 

The solution to this equation has to satisfy the first regularity condition in Eq. ([8), and the second 
boundary condition in Eq. (0. Using Eq. (fT3l l we rewrite the second regularity condition in terms 
of P, 



8P(S, ip) 



ds 



dP(s, -if) 



(16) 



s=o ds 

To derive the governing equations for non-axisymmetric tube oscillations we solve Eqs. (13[ 
and (TT3T > inside and outside the tube, and then match the two solutions at the tube boundary. It is 
straightforward to obtain the general solution to Eq. (fT3T > inside the tube satisfying the regularity 
conditions Eqs. © and ( TTol l. 

oo 

P' = ^ [cj, cosh(ns) cos(n<p) + D l n sinh(ns) sin(n^)] , (17) 

where O n and D), are arbitrary functions of t and z- The solution outside the tube has to decay as 
s — » oo. Hence, its general form is 

oo 

P e = J] e-" s [CI cos(n^) + D e n sm(n<p)] , (18) 

n=\ 

where once again C\ and D\ are arbitrary functions of t and z. Substituti ng Eqs. (TTTTt and (TT8l in 
Eq. (TT~3T > and using the expression for Vj_ in the elliptical coordinates (e.g. 



Korn & Kon 



^ = -=Fi\ 3 T + 9>^-)' = (sinh 2 S + sin 2 ^) 1/2 , (19) 
crfc) \ ds dip) 

where s and ip are the unit vectors in the s and ^-direction, we obtain the expressions for inside 
and outside the tube, 



1 °° 

£ = — — ) F) sinh(«i) cos(n^>) + G) cosh(ns) sin(ni,cj) , (20) 

O"0 L J 

n=\ 

1 °° 

^ = ~^ Z e ~" s [F " cos(w ^ + G " sin(n ^ )] ■ (2 1 } 



CT0 
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In these equations F' n , G l n , F e n and Gf, are functions of t and z- They are related to the functions O n 
D\, Q and 25* by 



u r n 2 u r n 






V . ■ = 

c\ f 2 Ai .3,2 
Or oz 


Pi 


(22) 


d G„ 2 d G„ 


J\ 




dfi Vai dz 2 


Pi ' 


(23) 


d 2 Ff, 2 #K 


c e 




dfi ' Ae 8z 2 


_ n 
Pe ' 


(24) 


d 2 Gl 2 d 2 Gl 


D e 




dfi lAe dz 2 


_ n 
Pe ' 


(25) 



Substituting Eqs. ( TTTl i and ( fl~8b in the second boundary condition in Eq. (O we obtain 

C n cosh(nso) = e-" Sa q v D), sinh(«s) = e-" s <>D e n . (26) 
Substituting Eqs. d20b and (1211) in the first boundary condition in Eq. © yields 

F\ sinh(ns ) = e _ " S(, F^ G\ cosh(ns) = e^G*. (27) 

Eliminating C'„, and F e n from Eqs. d22l) . d24l i. d26l and (f27T > we obtain the equation for FJ,, 

d 2 F n _ 2 d 2 /^ = 2 B 2 Q [1 + tanhQ» )] 
a<2 c »c &2 ' c nc w r j0i+ p e tanh(n*o)]' 

where we have dropped the superscript 'i'. Eliminating D\, D\ and GJ; from Eqs. ( l23b . (l25b . ( f26b 
and d27b we obtain the equation for GJ,, 

o 2 G 1 _ 2 ^ 2 gg[l+tanh(nj )] 

where we have once again dropped the superscript 'F. It follows from Eqs. Q and ( l20t that F„ and 
G„ have to satisfy the boundary conditions 

F n = 0, G„ = at z = +L/2. (30) 

In Eqs. d28l i and d29b « = 1 corresponds to kink modes, and n > 1 to fluting modes. 

In the elliptical coordinates the loop axis (x = y = 0) is defined by s = and ip = n/2. It 
follows from Eq. (1201 that the kink mode described by Eq. ( l28l i does not displace the loop axis 
in the ^-direction which, at the loop axis, coincides with the y-direction. Hence, the loop axis 
displacement is in the ^-direction, i.e. this mode is polarised in the direction of the larger axis 
of the tube cross-section. The kink mode described by Eq. ( f29T > displaces the loop axis in the s- 
direction. It is straightforward to show that it does not displace it in the ^-direction which, at the 
loop axis, coincides with the x-direction. Hence, the loop axis displacement is in the y-direction, 
i.e. this mode is polarised in the direction of the smaller axis of the tube cross-section. 

When the density is constant, we can use Eqs. ( 28b an d d29| i with the boundary conditions 



Ruderman 



(120031) . Let us look for the eigenmodes and 



Eq. (f30b to recover the results obtained by 
restrict the analysis to the fundamental modes in the z-direction. This implies that we take F„ and 
G„ proportional to e'" ' cos(nz/L). Then we immediately obtain that the eigenfrequencies of the 
boundary value problem defined by Eq. (l28T l and the boundary conditions ( l30l l are given by 

2 _ ^4 7r 2 ^[l+tanh(«so)] 
U ' K L 2 ^ 2 [p,+Petanh(«, )]' " ' """ ( ' 
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and the eigenfrequencies of the boundary value problem defined by Eq. d29b and the boundary 
conditions (l30t are given by 



2 



n 2 ci n 2 B 2 [l+tmh(n So )] 



n = 1,2, 



L 2 noL 2 [pi tanh(ns ) + p e ] 
In particular, the squares of eigenfrequencies of the kink modes are given by 
2 7T 2 Bl(a + b) 2 jr 2 Bl(a + b) 



/lc /i L 2 (api + bp e ) ' W ' s /ioL 2 (bpi + ap e ) ' 
It is straightforward to see that the eigenfrequencies satisfy 

LJ lc < a> 2c < ■ ■ ■ < Cl>2s < Wi s . 



(32) 



(33) 



(34) 



4. Implication on coronal seismology 



After 



Verwichte et al 



([2J304J) reported two cases of observations of the transverse coronal loop os- 



cillations where, in addition to the fundamental harmonic, the first overtone was also observed, 



Andries et al 



d2005l) sugg ested observations o f this nature could be used to estimate the scale 
(120051) assumed that an oscillating loop has a half-circle 



Andries et al 



height in the solar corona, 
shape and a circular cross-section, and it is in the vertical plane. They also assumed that the atmo- 
sphere is isothermal. In that case, the dependence of the plasma density on z is given by 

L 



p e = Pf exp 



/ L 7TZ\ 

cos — . 

\ nH LI 



^Pe 



(35) 



where H is the atmospheric scale height, p f the pl asma density at the loop foot points outside the 
loop, and £ > 1 a constant. 



Andries et al 



(120051) calculated the ratio of frequencies of the first 
overtone and fundamental mode and found that this ratio is a monotonically decreasing function of 
the parameter L/H. Hence, if we know the ratio of frequencies and L, we can determine H. F or a 
recent review of coronal seismology using kink oscillation over tones seelAndries et al.| (120091) . 
A very important que stion is how robust is this method. 



Dvmova & Ruderman 



(2006) and 



Morton & Erdelvil d2009al) have found tha t the account of th e loop shape can moderately affect the 
estimates of the atmospheric scale height. iRudermanl (120071) has shown that the twist of magnetic 
field lin es in the loop can be sa fely neglected when estimating the atmospheric scale height in the 



corona. 



Robertson et al 



(120101) found that the estimates of the atmospheric scale height obtained 



using the two-thread model are exactly the same as those obtained usin g the model of a m onolithic 



Ruderman 



(2010) showed 



coronal loop with a circular cross-section of constant radius. Recently 
that the account of stationary time independent siphon flows in coronal loops have little influence 
on the estimates of the cor onal scale height fo und using the frequency ratio. On the other hand, 



Ruderman et al 



(120081) and 



Verth et al 



(2008) found that the account of the loop expansion can 
strongly affect these estimates. 

In this section we study what the effect the elliptic cross-section has on the estimates of the 
coronal scale height. As we have already seen, when a loop has an elliptic cross-section, its kink 
oscillations are polarised along the axes of the cross-section. The kink mode polarised in the direc- 
tion of the larger axis is described by Eq. (l28T i with n = 1, while the kink mode polarised in the 
direction of the smaller axis is described by Eq. ( l29l l with n = 1 . Let us consider the solutions to 
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these equations in the form of eigenmodes and take F\ and Gi proportional to exp(-/w?)- Using 
Eq. d35l l we obtain 



B 2 Ja + b) 



Then, introducing 



at 



PoPf (at + b) 

2 = B 2 (a + b) 
ls ~ HQPfibt; + a) 

H<,Pi(a£ + b)a> 2 



exp 



B 2 Q (a + b) 

we reduce Eqs. d28l l and d29l with n = 1 to 

d 2 U 



exp 



— cos , 

\nH LI 

( 1 nz \ 

— cos 

\nH L ) 

_ poPiiH + a)co 2 
B 2 Q (a + b) 



arU 2 j L nz\ 
— — + QrU exp — cos — =0, 
dz 2 \nH L ) 



(36) 



(37) 



(38) 



(39) 



where either U — F\ and Q = f2 c , or U — G\ and Q, = Q s , and U satisfies the boundary conditions 
U — at z = +L/2. Since Eq. ( [391 does not contain a and b, the eigenvalues of the boundary value 
problem for U are independent of a and b. In particular, they are the same as those for a loop with 
the circular cross-section. Since 



£l 2 



^2s 



Q2 

Qi' 



it follows that we obtain the same estimates of the atmospheric scale height no matter if we use 
the observation of the kink oscillations polarised in the direction of the larger or smaller axis. The 
estimates are also independent of a and b and are the same as those obtained for a loop with the 
circular cross-section. 



5. Summary and conclusions 

In this paper we have studied non-axisymmetric oscillations of straight magnetic loops with a con- 
stant elliptic cross-section and density varying along the loop. We derived the governing equations 
for kink and fluting modes in the thin tube approximation. All these equations are similar to the 
equation describing kink oscillations of a straight tube with the circular cross-section. We found 
that there are two kink modes, one polarised in the direction of larger axis of the elliptic cross- 
section, and the other polarised in the direction of smaller axis. The frequencies of fundamental 
mode and overtones of these two kinds of kink oscillation are different. However, the ratio of fre- 
quencies of the first overtone and the fundamental mode is the same for both kink oscillations, and 
it is independent of the ratio of the ellipse half-axes a/b. This result implies that we obtain the same 
estimates of the atmospheric scale height no matter if we use the observation of the kink oscilla- 
tions polarised in the direction of larger or smaller axis. The estimates are also the same as those 
obtained for a loop with the circular cross-section. This demonstrates that the model shows a very 
robust nature when considering a static plasma. However, if the plasma in the loops is dynamic 
(i.e. time dependent) then the ability of the stati c model to provide accurate estimates may become 
questionable (see e.g. 



Morton & Erdelvi 2009b) 
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